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Abstract 

We consider the classical Wicksell corpuscle problem with spherical particles in M" and 
investigate the shapes of lower tails of distributions of 'sphere radii' in R" and 'sphere radii' 
in a fc-dimensional section plane Ek ■ We show in which way the domains of attraction are 
related to each other. 

1 Wicksell's corpuscle problem 

Suppose that a collection of spheres is randomly scattered in the Euclidean space M" and we 
intersect the collection with a fc-dimensional section plane E^, k G {1, . . . ,n — 1}. Then, under 
suitable model assumptions, the relation between the radii distribution F of spheres in and 
the radii distribution of corresponding spheres in the section plane is given by Wicksell's 

integral equation, 

= 1 _ r ^(^2 _ ^2^^(^l _ ^(^)) du, X>0. (1.1) 

where Mn-k £ (0, oo) denotes the (n — fc)th moment of F. Relation (|1.1|) has been discovered 
in the mid-twenties by S.D. Wicksell for the special case n = 3 and fc = 2. A proof of (jl.lj) 
can be found in Ch. 4.3]. Note that the density of F^"''^) always exists and is given by 

&^'\x) = ^-^^ Tiu'-xr-^dFiu), .>0, (1.2) 

Mn~k Jx 

which can easily be seen by differentiating (|1.H) . From (|l.lj) and (|1.2|) we see that the repre- 
sentations of and /("'^) only depend on the difference n — fc. In the following we will 
therefore write F^^'^ for f("'^\ /^^'^ for /("'^), and Mr for Mn-k, where r = n — k. 

In recent years there has been some effort in creating a theory of stereology of extremes. Drees 
and Reiss have investigated the shapes of upper tails of distributions of 'sphere radii' and 
'circle radii' that are connected by (|l.lj) . In this paper we investigate the lower tail behaviour. 

2 Extreme value theory and regular variation 

Let {Xi, i > 1} be a sequence of independent and identically distributed real valued random 
variables with common distribution function F. Denote the sample minimum by 

Wn = min(Xi, . . . , Xn) , n > 1 . 
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If there are sequences of normalizing constants a„ > and 6„ G M such that the normahzed 
minima — 6„)/fln converge in distribution to a nondegenerate distribution function H, then 
F lies in the domain of attraction of H. We denote this by -F G D^H). There are only three 
types of possible limiting distributions H, see e.g. |H|, namely 

|l -exp(-x") X > 0, i = 2, 

and 

iJ3(x) = 1 -exp(-e'^), xeR. 

Hi^ai -f^2,a and are called Prechet, Weibull and Gumbel distributions respectively. Recall 
the following condition for F G P(i?2,Q)- Let r/ = inf{x G M : F{x) > 0} denote the lower 
endpoint of F. Then F G T>{H2^a) if and only if > — oo and for all x > 0, 

lim^i;^=x°. (2.1) 

In the following we recall some basic facts from the theory of regular variation, see e.g. [H]. 
A function R is regularly varying at infinity with exponent p G M if it is real-valued, positive 
and measurable on [xq, oo), for some xq > 0, and if for each x > 

lim -z— — = x^ . 

t^oo R{t) 

We then write R G RVoo(/o). A function R{-) is regularly varying at if -R(^) is regularly 
varying at infinity with exponent (— /?) and this is denoted by G RVo(/o). If p = 0, is said 
to be slowly varying. A function R{-) is regularly varying if and only if it can be written in the 
form 

R{x) = xPL{x) , (2.2) 

where /? G M and L(-) is slowly varying. Note that the sum of two slowly varying functions is 
again slowly varying. For the investigation of the lower tail in Wicksell's corpuscle problem we 
need the following lemmas. 

Lemma 2.1. If Ri G RVo(/?i) and R2 G RVo(/32) then Ri + R2 £ RVo(min{/3i, ;g2}). 

Lemma 2.2 (cf. Theorem 2.7 in [5j). Let L be slowly varying at on (0, 00) and hounded on 
each finite subinterval of (0, 00). Suppose that for /? > the integral t^f{t) dt is well-defined 
for some given real function f and a given number 5 > 0. Then as u I 

00 />oo 

f{t)L{ut) dt ~ L{u) / f{t) dt 
for (5 > 0, and also for 5 = provided that L is non-increasing on (0, 00). 

Lemma 2.3 (cf. Exercise 1.1.1. in Let F G P(i?3) and r] = inf{x G M : F{x) > 0} > 
—00. Then for all n > 1, 

lim(x - riyF{x) = 0. 
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Lemma 2.4. Let /(•, •) and g{-, •) be positive functions such that 

j f{s, t) dt < oo , / g{s,t)dt <oo 
Jo Jo 

for some oj G (0, cxd]. Furthermore assume for s < t < uj 

lim "^^ ' } = c, c S [0, oo] . 
s^uj g[s,t) 

Then 

r f{s,t)dt 

Is9{s,t)dt 

3 Tail behaviour in Wicksell's corpuscle problem 

In recent years the main theoretical tools for estimating stereologically the tail of a particle size 
distribution have been stability properties of the domain of attraction. Consider the spherical 
Wicksell corpuscle problem with corresponding distribution functions F and F^'^\ r > 1, as 
described in Section ^ 

Note that there is no need to investigate the behaviour for F € T){Hi^a) since in this case 
r] = inf{x E M : F{x) > 0} = — oo. Let us first assume that r = 1 and F G P(ff2,a)- The radius 
of the spheres in M" is denoted by ^. We have 



Theorem 3.1. Let F £ P(i?2,a) and r] = 0. Then F^^) G V{H2,f3), where 



2, a > 1 , 

a + 1, 0<a<l. 



If ?? > 0, then e V{H2,2)- 

Proof. Let us first assume that r] = 0. It is easily seen from p.2j) that 

FW(t) = — / udF{u)+ {u - ^u^ - t^) dF{u) -- 
^1 [Jo Jt 

Using integration by parts in the first summand of (|3.1() we get 



:^[/i(t)+/2(t)]. (3.1) 



h{t) = tF{t)- f F{ 
Jo 



u) du . 



X 



as no, (3.2) 



For all a > and x > we show that 

h{tx) 
hit) 

i.e. Ii G RVo(a + 1). For that consider 

hitx) _ txF{tx) - /o"^ F{u) du _ xF{tx) 1 - {txF{tx))-^ /q"' F{u) du 
hit) ~ tFit) - /o Fiu) du 1-itFit))-^ J^Fiu)du 

Define G(s) = s"^F(i) G RVoo(— a — 2) and apply Karamata's theorem Th. 0.6] to obtain 

^^ f^Fiu)du rGis)ds 1 



no tFit) z^oo zGiz) a + 1 ' 



3 



which proves ()3.2() . 

For the second summand l2{t) we consider three cases. 
Case I, a > 1: Since a > 1, E(^~-'^) < oo. Rewriting /2(i) as 



hit) = E ((e - ye^)%>t}) = t'B ((e + v'e^)-'%M}, 

and applying the dominated convergence theorem we get 

hmE ((e + Ve-t^r'l{i>t}) = ^E(r') G (O, oo) , (3.3) 

which imphes h G RVo(2). By Lemma O F^^) € RVo(2), whence F^^) G V{H2,2)- 

Case II, a G (0, 1): We show that I2 G RVo(a + 1), which imphes F^^) G RVo(a + 1) by Lemma 

12.11 For that we use again integration by parts to arrive at 

hit) = -tF{t) - / F(u)(l - ^-^ ) dn . (3.4) 

Writing F(t) = t'^L[t) where L is slowly varying at 0, and substituting u = tv \n the integral of 
((231) yields that 

/oo 
L(ti;)w°g(t;) dt- , 

where 

g{v) = ^^-l, v>l, (3.5) 

is a probability density function. For (5 G (0, 1 — q) we have th'^g{t) dt < 00 and therefore 
Lemma 12.21 implies 

f?^ L(tv)v°'q(v) dv 

Ht):=- ^ ' ^1 astiO, 

cL(t) 

where c := v'^g{v) dv G (1, 00). Thus 

h{t) = t''+^L{t)[ch{t) - 1] G RVo(a + 1) . 

Case III, a = 1: We show that h G RVo(2). Choose A > 1 such that c{A) := jf vg{v) dv > 1. 
From (|3.4|) we get 



hit) = t 



F{vt)g{v)dv - F{t) + F{vt)g{v)dv 



A 



(3.6) 



where g is given by (|3.5jl . Consider 



/A pA 
F{vt)g{v) dv - F{t) = t J L{vt)vg{v) dv - tL{t) , 

where L is slowly varying at 0. Using the uniform convergence theorem for slowly varying 
functions we get 

, * / N L{tv)vg{v) dv 

h*{t):= ^^ ^1 as no, 

c{A)L{t) 

whence 

Ji(t) = tL{t)[c{A)h* (t) - 1] G RVo(l) . 
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The second summand in ()3.6p can be written as 

/"OO /"OO 

J2{t) = / F{vt)g{v) dv = / V{v-^i < t)g{v) dv = P(r'e < t, C > ^) , 
J A J A 

where C is a random variable with probabihty density function g independent of ^. Furthermore 

P(C- < t, C > ^) = ^{CH <t\C> A)V{Q >A) = F~^,,^{t)V{C > A) , 

where C is a random variable with probability density function g : {A, oo) — > (0, oo) given by 
g{u) = c~^g{u) for c = g{u) du. It is easy to check that G RVo(l). By assumption 

G RVo(l) and since ^ and C,~^ are independent we have -^^-i^ £ RVo(l) (see e.g. Theorem 
3 and Corollary therein) and therefore J2 G RVo(l). Hence 

h{t) = t[Ji{t) + j2{t)]^mo{2). 

Now let r/ > 0. For sufficiently small t we get from 

For all a > we have E(^^^) < 00, whence F^^^ G 'D{H2^2), using the same argument as in 
(|3.Hj) . The proof is complete. □ 

The case of a Gumbel limiting distribution is considered in the next theorem. 
Theorem 3.2. If F £ V{H-i), then F^ g V{H2,2)- 

Proof. The case r/ > is covered by Theorem 13.11 So let rj = 0. From 1)3. ip we see that 

where 771 is a random variable with density g given by (|3.5j) and r/2 is uniformly distributed on 
[0, 1] and 7/1 , , are independent. It suffices to show that -F^/^^ G RVo(l) and F^/j^^{t) / F^/^^{t) — 
as t i 0. 

Let us first prove that F^/^^ G RVo(l). For that we write 

F^/^^it) = P(e < tm,vi < lA) + P(^ < tvuvi > lA) 

1 1 



F{u)g{u/t)du + - I F{u)g{u/t)du 
t Jt t Ji 

=: J{{t) + r^{t) . (3.7) 
For the second summand in (|3.7|) we can apply Lemma 12.41 and get 

j^{tx) lyt ntw){^^^ - 1) dw 

Since g G RVoo(— 2) the first summand in (|3.7|) can be written as 

J*{t)=t[ F{u)u'^L{u/t)du, 
Jo 



5 



where L{l/t) is slowly varying at 0. 



Because of Lemma \2.'A\ u ^ '^F(u) du < oo for some 6 > and therefore we can apply [HI 
Theorem 2.7] and get 



h**(t) = ^ — ^—^ .1 astjO, 

L{l/t)J^Fiu)u--^du 



J^{t) = tL{l/t)h**{t) I F{u)u-^ du G RVo(l) . 
Jo 

Next we show that F^/jj,^{t)/ F^/j^^{t) ^ as t | 0. Since G RVo(— 1), we have 

limi~^Ft/„, (t) = oo . 
no ^' 

Furthermore Lemma 12.31 implies 

limt~^Ff:/„Jt) = limt"^ [ F(tu)du = 0. 
HO '^"^^ ' no Jo 

The result follows. □ 

The sectioning of the system of spheres with a /c-dimensional plane can be obtained by an 
iterated intersection procedure. First we intersect the system with an (n — l)-dimensional plane 
En-i, then the obtained system of spheres in En-i is intersected with an (n — 2)-dimensional 
plane and so on. This leads to 

Corollary 3.3. Let F G P(i?2,a). Then G P(i?2,2) for r > 2. If F e V{H-i), then 
Fir) g x>{H2,2) for aii r > 1. ' 
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